Resistance Genes, Plasmids, Multilocus Sequence Typing (MLST), and Phenotypic Resistance of Non-Typhoidal Salmonella (NTS) Isolated from Slaughtered Chickens in Burkina Faso

The emergence of antimicrobial-resistant bacteria in developing countries increases risks to the health of both such countries’ residents and the global community due to international travel. It is consequently necessary to investigate antimicrobial-resistant pathogens in countries such as Burkina Faso, where surveillance data are not available. To study the epidemiology of antibiotic resistance in Salmonella, 102 Salmonella strains isolated from slaughtered chickens were subjected to whole-genome sequencing (WGS) to obtain information on antimicrobial resistance (AMR) genes and other genetic factors. Twenty-two different serotypes were identified using WGS, the most prevalent of which were Hato (28/102, 27.5%) and Derby (23/102, 22.5%). All strains analyzed possessed at least one and up to nine AMR genes, with the most prevalent being the non-functional aac(6′)-Iaa gene, followed by aph(6)-Id. Multi-drug resistance was found genotypically in 36.2% of the isolates for different classes of antibiotics, such as fosfomycin and β-lactams, among others. Plasmids were identified in 43.1% of isolates (44/102), and 25 plasmids were confirmed to carry AMR genes. The results show that chicken can be considered as a reservoir of antibiotic-resistant Salmonella strains. Due to the prevalence of these drug-resistant pathogens and the potential for foodborne illnesses, poultry processing and cooking should be performed with attention to prescribed safe handling methods to avoid cross-contamination with chicken products.


Introduction
Foodborne illness caused by Salmonella is a public health concern around the world [1]. Salmonella are zoonotic bacteria and show increasing resistance to antibiotics, especially those isolated from farm animals. Several researchers have shown that Salmonella circulate between animals, humans, and the environment [2][3][4]. Chickens constitute a reservoir of Salmonella, which can be part of the normal intestinal flora of these birds [5]. In Burkina Faso, chicken farming is a booming industry, and the consumer demand for chicken is increasing every year. In Ouagadougou (population~3 million), the capital city of Burkina Faso, it is estimated that more than 80,000 chickens are slaughtered every day for consumption [6]. Chicken farming has become a good source of income for many Burkinabe who are carrying out this activity without adequate training on good farming hygienic practices. In this country, animal feed is produced in an unregulated manner without veterinary checks or microbiological quality analysis. In addition, there are feed producers who add antibiotics to their feed to prevent disease and the loss of flocks. These practices could select for antibiotic resistance in the bacteria of the chickens' intestinal flora. This has greatly affected the emergence of multi-drug-resistant strains of Salmonella. The slaughter conditions of chickens constitute a critical point of contamination of carcasses, especially during evisceration [7]. After slaughter, the carcasses are exposed to ambient temperature for sale all day. These practices undoubtedly constitute the critical points of the contamination of humans by pathogenic bacteria from chickens [7,8]. Previous data have shown the presence of antibiotic-resistant Salmonella in Burkina Faso, but very little data exist on the genomic characteristics of these Salmonella due to a lack of equipment and methods in Burkina Faso [9,10]. Whole-genome sequencing (WGS) allows for the analysis of the Salmonella genome but requires expensive equipment and reagents, which are inaccessible to laboratories in Burkina Faso; therefore, for this study, we collaborated with scientists at the United States Department of Agriculture, Agricultural Research Service (USDA-ARS). Previously, we investigated the genomic characterization of Salmonella isolated from fish in Burkina Faso [11]. The present study aims to characterize Salmonella strains isolated from slaughtered chickens in Burkina Faso using WGS to better understand their molecular epidemiology.

Serotypes of Isolates from Slaughtered Chickens
Twenty-two different serotypes were identified using WGS data. The most prominent serotype found was Hato with 28 (27.5%) isolates, followed by Derby with 23 (22.5%); Muenster with 7 (6.9%); and Typhimurium, Poona, Chester, and Kentucky all with 4 (3.9%). Other serotypes found include Alexanderplatz and Bredeney with three isolates (2.9%) each, and Rechovot, Telelkebir, and Tennessee with two (2.0%) each. Five isolates (4.9%) were named with their antigenic formula, and one isolate serotype could not be distinguished between Albany and/or Dusseldorf (Table 1 and Figure 1).
Multiple resistance genes were found in 37 (36.3%) of the strains, with each strain possessing two to nine genes. Nineteen (19/28) S. Hato isolates possessed seven to nine resistance genes, including up to five aminoglycoside resistance genes, sul2 conferring resistance to sulfonamides, and tet(A) conferring resistance to tetracycline, and two contained dfrA14 conferring resistance to trimethoprim. Twenty-one of the S. Hato isolates possessed the sul2 gene, and nineteen had the dfrA gene. Twenty of these isolates showed phenotypic resistance to trimethoprim/sulfamethoxazole. Of the 19 isolates possessing the tet(A) gene, 18 showed phenotypic tetracycline resistance. However, in 15 of the 28 S. Hato isolates, at least one AMR gene was a partial sequence (Table 1).

Antimicrobial Resistance Genes and Antibiotic Resistance Phenotypes Detected
Antibiotic resistance genes and phenotypic resistance profiles are shown in Table 1 and Figure 1. Phenotypic resistance to at least one antibiotic was observed among 33 Of the 23 S. Derby isolates, 7 contained multiple resistance genes. Five isolates had five resistance genes, and one isolate had six resistance genes, including multiple aminoglycoside genes, sul2, tet(A), and fosA7. Two S. Derby isolates had the fosA7 gene. The six S. Derby isolates with tet(A) were phenotypically resistant to tetracycline, with four also being resistant to Minocycline.
All seven S. Muenster isolates had only one resistance gene, aac(6 )-Iaa, but no aminoglycoside resistance was seen. Of the four S. Typhimurium isolates, one possessed bla TEM-1B , conferring resistance to β-lactams, and mph(A), suggesting resistance to macrolides. However, no macrolides were tested on this panel. Interestingly, one S. Typhimurium isolate with only the aac(6 )-Iaa gene found by WGS showed resistance to tetracycline and trimethoprim/sulfamethoxazole. Of the four S. Kentucky isolates, three had nine resistance genes conferring resistance to aminoglycosides, sulfonamides, trimethoprim, tetracyclines, and quaternary ammonium compounds (qacE). The qacE gene was partially present in all S. Kentucky strains, which may suggest resistance to sanitizers such as benzalkonium chloride. The chromosomal mutations and MLST results of the strains are shown in Table 1. Twenty-eight different MLST types were identified. Multiple MLST types were identified among serotypes Derby, Hato, Poona, Telelkebir, and Typhimurium. Some isolates of the serotypes Alexanderplatz, Brancaster, Derby, Farmingdale, Hato, I 1,3,19:b:-, I 1,3,19:f,g:1,5, and Rechovot had unknown MLST types.
Among the isolates with multiple AR genes, strong positive correlations (r > 0.7) were found for the co-occurrence of several AR genes ( Figure 2). Strong positive correlations were found for the co-occurrence of multiple aminoglycoside genes, including aadA1 with aph(3 )-Ia, aph(3 )Ib, or aph(6 )-Id. Strong correlations were also found for the co-occurrence of antibiotic resistance genes from different antibiotic classes, including sul2 with multiple aminoglycoside resistance genes, tetA, or drf A14 ( Figure 2).
Plasmid replicons and AMR genes were present on the same contig in 25 isolates (Table 2). However, other AMR genes were likely physically linked to plasmid replicons but were not detected in the WGS assembly as evidenced by some isolates containing AR genes and plasmid replicons with strong positive correlations (r > 0.7) for the cooccurrence of AR genes and plasmid replicons. In particular, the IncQ1 replicon was strongly correlated with tetA, sul2, dfrA14, and multiple aminoglycoside resistance genes ( Figure 2 and Supplementary Table S1).

Phylogenetic Analysis of S. Derby Isolates
Both S. Derby isolates containing the fosA7 gene were compared to publicly available S. Derby genomes from chicken sources using cgMLST (Figure 3). The two isolates from this study were the only members of cgMLST type 227,637 in the analyzed dataset and were located on a branch by themselves.

Discussion
The present study shows that poultry is a reservoir of MDR Salmonella strains in Burkina Faso. Kagambèga et al. [4] previously reported this. A previous study demonstrated that Burkina Faso does not have a commercial slaughterhouse for chickens and that slaughtering is instead carried out at traditional markets [9]. The conditions of chicken slaughter do not respect good hygienic practices, and this undoubtedly promotes cross-contamination, especially during evisceration, between carcasses and chicken feces, representing a risk to human health in Burkina Faso. Figure 3. Phylogenetic analysis of the examined S. Derby isolates containing fosA7 (n = 2) and publicly available S. Derby isolates from chicken in Enterobase (https://enterobase.warwick.ac.uk/ accessed on 1 June 2022) using single-nucleotide polymorphisms (SNPs) and hierarchical clustering of core genome (cg) MLST (HierCC). The legend shows cgMLST HC100, which indicates allelic differences of no more than 100 of 2850 core genomic alleles among isolates.

Discussion
The present study shows that poultry is a reservoir of MDR Salmonella strains in Burkina Faso. Kagambèga et al. [4] previously reported this. A previous study demonstrated that Burkina Faso does not have a commercial slaughterhouse for chickens and that slaughtering is instead carried out at traditional markets [9]. The conditions of chicken slaughter do not respect good hygienic practices, and this undoubtedly promotes cross-contamination, especially during evisceration, between carcasses and chicken feces, representing a risk to human health in Burkina Faso.
MDR Salmonella Typhimurium and Kentucky regularly cause human salmonellosis in Burkina Faso, but the sources of these salmonellosis cases remain uninvestigated [13,14]. Many researchers have demonstrated that poultry eggs and meat are major vehicles for human salmonellosis, which is exacerbated by imports from around the world [4,15,16]. The uncontrolled use of antibiotics in Burkina Faso poultry farming contributes to the development of AMR in pathogens and commensal flora. Moreover, the emergence of AMR in bacteria from poultry farms has generated human health concerns due to the consumption of contaminated meat and eggs [17].
Salmonella Hato and Derby were the most prevalent serotypes isolated from chickens in this study. Previous studies in Burkina Faso revealed that S. Derby and S. Hato were the major serotypes circulating in poultry [9,10]. Abdelkader et al. [18] found similar results in Niger, which borders Burkina Faso. S. Hato and S. Derby isolates showed phenotypic resistance to different classes of antibiotics, such as aminoglycosides, tetracycline, trimethoprim, and sulfonamides, and contained genes predicted to confer resistance to fosfomycin. Interestingly, eight S. Derby and two S. Hato isolates contained a phage-like plasmid with 99% identity to AnCo3. Phage-like plasmids have previously been described in North America to be associated with bla CTX-M-15 genes. AnCo3 was first identified in a S. Derby clinical isolate in Canada. Although similar to AnCo and AnCo2, which contain bla CTX-M-15 , this phage-like plasmid contains no AMR genes. To the best of our knowledge, this is the first report of a phage-like plasmid in Burkina Faso, indicating the global spread of these emerging mobile genetic elements.
However, several of the AMR genes identified in S. Hato isolates were only partial sequences. Despite the genes for sulfonamide and trimethoprim resistance appearing as partial in the WGS data, most of these isolates still showed phenotypic resistance to trimethoprim/sulfamethoxazole. These isolates also contained partial IncQ1 replicon sequences, some of which were co-located on contigs with IncI1 replicons. It is possible that these plasmids have merged and that the AMR genes carried have been disrupted. It is also possible that these isolates contained yet unknown resistance genes for these antibiotics. However, it is also possible that an assembly error stemming from repeated DNA sequences, which is notoriously difficult to assemble, caused these genes to appear partial when they are in fact not. In all of these cases of partial sequences, long-read sequencing would be beneficial to further investigate the genetic structure of these plasmids and AMR genes.
Unsurprisingly, the S. Derby isolates containing the fosA7 gene in this study were genetically unique as compared to publicly available sequences from chickens. The majority of the publicly available Salmonella genomes are from countries with surveillance systems, which Burkina Faso lacks, so it is reasonable that these geographically distinct isolates would also be genetically unique.
Salmonella Typhimurium isolated in this study possessed bla TEM-1B β-lactamase for β-lactam resistance. Extended-spectrum β-lactam resistance could not be confirmed on the susceptibility panel. One S. Typhimurium isolate possessed the resistance gene mph(A) for macrolide resistance, which could not be confirmed because the panel lacked a macrolide antibiotic. MDR S. Typhimurium strains were previously isolated from chickens, and they showed more than 80% genetic similarity to S. Typhimurium isolated from human patients [4]. These facts show that chickens and their products constitute a potential danger for the colonization of humans with antibiotic-resistant Salmonella. In this study, 36.3% (37/102) of the strains analyzed contained resistance genes for two or more antibiotic classes, with resistance genes from the aminoglycoside class being the most prevalent. Some strains contained up to nine different resistance genes.
The resistance genes found in this study did not always correlate with the resistance phenotypes. While many isolates possess aminoglycoside resistance genes, the antimicrobial susceptibility test results showed that only four (3.9%) isolates were resistant to gentamicin and that none were resistant to tobramycin. One S. Typhimurium isolate only possessed the aac(6 )-Iaa gene but showed phenotypic resistance to tetracycline and trimethoprim/sulfamethoxazole. In this case, it is possible that the strain lost a plasmid containing the genes for resistance to these drugs between susceptibility testing and sequencing or that the isolate contains yet unknown genes.
Several antibiotic resistance genes were also detected where phenotypic resistance was not confirmed. For example, there is no CLSI method for phenotypic determination of fosfomycin resistance using broth microdilution, so phenotypic resistance could not be determined. Additionally, AAC(6 ) enzymes inactivate aminoglycoside antibiotics by acetylating their substrates at the 6 position and can confer resistance to amikacin and kanamycin, which were not included on the panel of antibiotics used [19]. However, the aac(6 )-Iaa gene has been demonstrated to be non-functional in Salmonella unless the strain possesses a mutation to render the promotor for the gene functional [20].
Three of the four S. Kentucky isolates identified in poultry in this study showed multidrug resistance phenotypically, with nine different resistance genes conferring resistance to four or more different antibiotic classes and the qacE gene conferring resistance to antiseptics, although this gene was only partially present. In contrast, a study conducted by Chuanchuen et al. [21] in Thailand found that twenty-seven percent of the Salmonella strains isolated from poultry and swine possessed qacE∆1 and that none of them harbored qacE. This qacE identified in this study could be explained by the repeated usage of disinfectants, including quaternary ammonium compounds (QACs), in the farm environment in Burkina Faso. This may increase the selection and persistence of bacteria with reduced susceptibility not only to antiseptics but also possibly to antibiotics [22]. However, phenotypic resistance would have to be confirmed as only a partial gene is present.
The use of PlasmidFinder in this study detected the plasmids with replicon sequences IncFIB(S), IncFII(S), IncFII(pCoo), IncFIB (H89-PhagePlasmid), IncB/O/K/Z, IncX1, and IncQ1 in Salmonella isolates. These isolates carried resistances genes for four, three, and two/one classes of antibiotics, including aminoglycosides, β-lactams, sulfonamides, tetracyclines, and phenicols. Villa et al. [23] reported similar results on the IncF group carrying ESBL or plasmid-mediated quinolone or aminoglycoside resistance genes. Carattoli et al. [24] demonstrated that the IncF plasmid family is prevalent in clinically resistant isolates of Enterobacteriaceae. Moreover, IncF can be virulence-associated plasmids, which give host bacteria the ability to cause a more virulent infection [25]. IncX1 plasmids have been associated with genes for resistance to β-lactams and aminoglycosides in Salmonella isolated in the USA and genes for resistance to quinolones globally [26,27]. Isolates carrying IncQ1 replicons also carried aminoglycoside resistance genes and sul2 with a strong correlation. This result was not surprising, as IncQ1 plasmids are known to be commonly associated with genes for resistance to aminoglycosides, tetracyclines, and sulfonamides [28].
The present study concurs with previous research that the T57S substitution detected in parC is not always associated with a quinolone resistance phenotype since it has been found in both resistant and susceptible isolates [11,29]. Feng et al. [30] also found similar results in a study of a Salmonella Goldcoast lineage in Northern Taiwan, where a single T57S mutation was not always sufficient to confer clinically significant resistance.

Bacterial Strains
Isolates were collected during a previous investigation of Salmonella found in various foods, food animals, and humans in Burkina Faso [9]. The Salmonella isolates (n = 102) from the cecal and/or intestinal contents of slaughtered chickens used in this study were obtained from the Laboratoire de Biologie Moléculaire, d'épidémiologie et de surveillance des bactéries et virus transmissible par les aliments (LaBESTA)/Université Joseph KI-ZERBO, Burkina Faso. Slaughtered chickens were sourced from markets in different villages across the country. Salmonella were isolated using standard methods as previously described [9].
4.3. DNA Extraction, Whole-Genome Sequencing, Assembly, Annotation, and Molecular Serotyping DNA extraction, library preparation, whole-genome sequencing, assembly, and annotation for the 102 Salmonella strains were completed as previously reported [11]. Briefly, libraries were prepared using Nextera XT DNA library preparation kits, which were sequenced using either a 300 or 500 cycle Illumina MiSeq version 2 reagent kit. Reads were assembled using A5 and annotated with the NCBI Prokaryotic Genome Annotation Pipeline [32]. The sequences were deposited into NCBI under BioProject no. PRJNA679582 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA679582 accessed on 22 March 2022). The serovar determination of the strains using SeqSero was previously described [11]. Identification of serotypes, antibiotic resistance genes, chromosomal mutations, MLST, and plasmids was carried out.
Antibiotic resistance genes and chromosomal point mutations associated with resistance were identified using ResFinder 4.1 through the Center for Genomics Epidemiology (CGE) website (https://cge.cbs.dtu.dk/services/ accessed on 6 May 2022) [33]. Genes with 80% identity and greater than 90% coverage were considered present, and genes with coverage between 40% and 90% were considered present but partial genes. Complete but disrupted genes were noted. Multilocus sequence type (MLST) for each isolate was identified using MLST 2.0 through CGE [34]. PlasmidFinder 2.1 accessed through CGE was used to identify plasmid replicons [35,36]. Replicons with 80% identity and greater than 90% coverage were considered present, and replicons with coverage between 40% and 90% were considered present but partial replicons. Partial AMR genes and replicons were confirmed as partial or complete but disrupted using BLAST [37]. For incompatibility groups with an established scheme, plasmid MLST (pMLST) type was determined using https://pubmlst.org/organisms/plasmid-mlst/ (accessed on 6 May 2022) [38]. Plasmid replicons or pMLST gene targets found on the same contig as AMR genes were noted. Contigs containing phage-like plasmid replicons were confirmed as phage-like plasmids with BLAST.

Phylogenetic Analysis of Salmonella Derby Isolates
Raw paired-end fastq files of both S. Derby isolates (S171 and S251), which contained fosA7, were imported into Enterobase (https://enterobase.warwick.ac.uk/, accessed on 6 May 2022) and compared to all the publicly available genomes of S. Derby (n = 197) sourced from poultry in Enterobase, updated on 12 October 2021, using single-nucleotide polymorphisms (SNPs) and hierarchical clustering of core genome (cg) MLST (HierCC) (Zhou et al., 2020). Our study isolates (S171 and S251) and the retrieved genomes from Enterobase were all aligned to the reference S. Derby 2014LSAL01779 complete genome (CP026609.1) and designated to HC100 differing by ≤100 core genomic alleles.

Statistical Analysis
To determine the overall distribution of plasmid replicon types and antimicrobial resistance genes among the examined Salmonella serovars, a heatmap with hierarchical clustering was generated using package "pheatmap" in R software (version 3.4.2). A correlation analysis was also performed to determine the association of both determinants among the examined Salmonella isolates. Antimicrobial resistance genes and plasmid replicons results, including partial sequences, were converted into binary data (0/1), where the presence of plasmid replicons and resistance genes in isolates received scores of 1, whereas absence of both determinants received scores of 0. The binary data (0/1) for antimicrobial resistance genes and plasmid replicon types were uploaded into R software (version 3.6.1; https://www.r-project.org, accessed on 6 May 2022), and the correlation was calculated at a significance of p < 0·05 using "cor" and "cor.mtest" functions. The correlation plot was then generated using the "corrplot" function. Based on the values of r, the degree of correlation is considered strong, moderate, and weak if r value is >0.6, 0.4-0.6, and <0.4, respectively.

Conclusions
This study shows once again that chicken constitutes a reservoir not only of pathogenic bacteria but also of bacteria that are multi-drug resistant. Chicken is a good source of animal protein and is very popular in Burkina Faso, as production is expanding in the country. This food-producing animal is a reservoir of multi-drug-resistant Salmonella, but this study is the first one in the country that reports potential resistance to antiseptics associated with multi-drug resistance. Unfortunately, the country does not have suitable slaughterhouses for chickens, and each market designates one site to slaughter, sell, and roast chickens, which very often smells foul and is visibly unsanitary. The soil and detritus from these poultry processing sites are transported by rainwater to the environment, which contributes to the pollution of water reservoirs and the environment with resistant pathogenic bacteria. Authorities should consider implementing a farm-to-fork quality control system to minimize the risk of pathogenic bacteria contaminating chickens and, subsequently, consumers. Using WGS is a very quick solution to characterize the genome of pathogenic bacteria. Efforts must be made to popularize these methods in developing countries such as Burkina Faso.